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Abstract. We study a singular parabolic equation of the total variation type in one dimension. The 
_C ■ problem is a simplification of the singular curvature flow. We show existence and uniqueness of weak 

solutions. We also prove existence of weak solutions to the semi-discretization of the problem as 
well as convergence of the approximating sequences. The semi-discretization shows that facets must 
form. For a class of initial data we are able to study in details the facet formation and interactions and 
their asymptotic behavior. We notice that our qualitative results may be interpreted with the help of a 
k>" ' special composition of multivalued operators. 
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1 Introduction 



Many free boundary problems involving the Gibbs-Thomson relation may be considered as a driven 
weighted mean curvature flow coupled through the forcing term to a diffusion equation (see MCRI1 , 
[Ra], [L] [AW]). We have a considerable body of literature concerning this problem for the Euclidean 
curvature of the interface, including the question of precise regularity of solutions treated by Escher, 
Priiss, Simonett and Mucha, see [EPSl, lESL |Mu]. On the other hand, less is known if the curvature 
appearing in the Gibbs-Thomson relation is singular, see e.g. |Ry[ . This line of research has been 
initiated by Taylor, [T], and independently by Gurtin, MAGI . However, just solvability of equations of 
the singular curvature flow is interesting. Existence of the flow was obtained by Bellettini, Novaga, 
Paolini [BNP1], [BNP2] and by Chambolle BChl . Driven singular curvature flow was studied by 
M.-H.Giga, Y.Giga and Rybka, see flGGl, llGRTTl . llGR2l . 

In fact, the existence and properties of solutions to the singular weighted mean curvature flow 

V = k on T(t), (1.1) 

are interesting in itself even in the plane and without forcing, especially when the anisotropy function 
(also called 'energy density function') is singular, i.e. just convex. Here, T(t) is the unknown curve 
and k denotes the weighted mean curvature related to the underlying anisotropy function and V is the 
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velocity of surface T(t). Our ultimate goal would be to study existence and behavior of solutions to 

CCD. 

In its full generality problem dl.ll ) for an arbitrary initial curve is rather difficult. One source of 
difficulties is the geometry of the system, it is already present in the two-dimensional setting. Here, 
we want to concentrate only on the purely analytical difficulties appearing in dl.ll ). This is why we 
will restrict our attention to a simplified equation, which retains the singular character of the original 
problem. 

Here is our postulated equation 

A t = |--^J( S + A s ) in,Sx(0,r), 
as a<p 

A(s,0)=A (s) onS", (1.2) 
A(2vr,t) = A(0,t), t>0, 

here S is the unit circle parameterized by interval [0, 2tt) and A is the sought function. Compared 
with dl.ll) our new system has one analytical advantage. Namely, the domain of definition of A(-, t) 
is independent of time. 

We present a justification of this equation in the Appendix. Here, we explain our notation. The 
variable s plays the role of the arclength parameter, the subscript s denotes the differentiation with 
respect to s. We frequently refer to <p = A s + s as the angle between the x\ axis and the outer normal 
to the curve. Such an interpretation helps drawing pictures, but the relation to the actual angle is rather 
loose. 

We make a specific choice of J corresponding to the surface energy density functions. We want 
to study a situation which is already very singular yet tractable. In many instances of a great physical 
interest an anisotropy appears, which is merely convex, not even strictly convex (understood in a 
proper sense). As a result, we choose J, which is convex and piecewise linear. This is an independent 
source of difficulties. In order to avoid further technical troubles we will choose J corresponding 
to the situation where that curve minimizing the surface energy (which is the Wulff shape of the 
anisotropy function) is a square. We must stress again that the correspondence is at the level of ideas, 
because (1 1 -2b is not a curvature flow, but its caricature. However, the obtained behavior of solutions 
to (1 1 -2b is almost the same as for the equation dl.ll ) with the anisotropy function corresponding to a 
square, HChl . 

Thus, we pick J which suffers jumps of equal height | at the equi-spaced angles 

A = la k = — — + kAa: k = 0,1,2,3, with Aa = - \. (1.3) 
Specifically, we put 

m = 4 (\<P ~ T l + \<P ~ 4 I + W + 4 I + W + T \) ■ (1-4) 

Since A is defined over the unit circle its graph over S is a closed curve. The meaning of the spacing 
between a^'s can be explained by looking at the equation 

--J(, + A s) = l, 

considered in | |MRy[ - see subsection 3.2, too. Roughly speaking, the spacing between a& and a^+i 
corresponds to the length of facets having the normal vector n with the normal angle a^. The size 
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of the jump of £j(s + A s ) corresponds to the angle between the normals to the curve, which is a 
solution to the above equation, at a corner. 

The chosen anisotropy function (11.41) is nowhere regular, hence we can expect nonstandard effects 
requiring new analytical tools. This has been observed by researchers working on the total variation 
flow, whose simplification is 

u t - S (u x )u xx = (1.5) 

augmented with initial and boundary data. Here, 5 a is the Dirac measure concentrated at a. 
We noticed so far two main types of motivation to study the total variation flow, 

u t- div {w\) =0 - (L6) 

The first one is the image denoising and reconstruction introduced by Rudin and Osher, [RO], [ROF]. 
The second one is evolution of the facets of crystals. The bulk of the papers (see IABC1I . [ABC2], 
[BCN], HGK1 . BGGKL BACM1 . HMoIO uses the theory of nonlinear semigroups to establish existence. 
The last paper is particularly interesting, because it deals with the anisotropic total variation flow. 
Moreover, the notion of entropy solutions was introduced to deal with uniqueness of the total variation 
flow (see MABC1B . IIBCNI ). The tools of convex analysis were useful to make sense out of (1 1 -5b - 
The authors, mentioned above, paid special attention to piecewise constant initial data and they were 
interested in the asymptotic behavior, in particular the asymptotic shape was identified. M.-H.Giga, 
Y.Giga and R. Kobayashi, MGKl . IIGGKB . also calculated the speed of flat facets. No matter what is 
the approach, it is apparent that the most important information is located in sets {u x = 0}, where the 
singular dissipation starts to play a role and where the classical multivalued theory of function loses 
the meaning. 

Our approach differs in many aspects. We prove existence by a regularizing procedure and passing 
to the limit with the regularizing parameter, this approach was used, e.g. by Feng and Prohl, see HFPI1 . 
The main difficulty is associated with studying the limit of the non-linear terms. We present a more 
detailed analysis of regularity of solutions permitting us to call them 'almost classical'. For generic 
data, our solutions are twice differentiable with respect to s, except a finite number of points (for fixed 
time). This will be explained in detail below. We mention here that we use the tools of the convex 
analysis, in particular we rely on the fact that for a convex function the subdifferential is well-defined 
everywhere. However, the classical theory of multivalued functions is not sufficient. We have to 
introduce a new definition of the composition of two multivalued functions to describe the meaning 
and qualitative properties of solutions to system (11.2b as well a class of the J-R functions, where 
regularity is described from the point of view of the properties of the function J. In our opinion the 
results we prove contribute to better understanding parabolic systems with measure coefficients. 

Our technique requires a new look at the regularity of functions. We will generalize the meaning 
of the convexity defining a class of J-regular functions preserving some important properties of the 
convexity. Our main qualitative result says that any sufficiently regular initial curve evolving accord- 
ing to system (11.21 ). will eventually reach a minimal solution, which is called the asymptotic profile 
in the area of the total variation flow. The geometric interpretation is that the solution reaches its 
asymptotic shape, i.e. the square in our case. This may happen in infinite of finite time depending 
upon initial data. If this event occurs in finite time, then subsequently, the solution shrinks to a point. 
This behavior can be illustrated by the pictures below. The precise meaning is contained in Theorem 

ED 
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Time t = Time t = t\ 

The evolution is determined by motion of facets defined by singularities of the J-function (the 
arrows show the direction of the evolution). In finite time we obtain a convex domain, which becomes 
a square converging to a point in finite time. 




Time t = t% Time t = t% Time t = t en d 

We have to underline that the illustrated evolution hides the novel idea of definition of singular 
term 5q{u x )u xx being a multiplication of two Dirac deltas (as in (11.5I )). however the nonlocal character 
will allow us to define this object. Additionally, by the uniqueness of solutions to our system we show 
that our novel definition is the only admissible. Formally, the dissipation caused the Dirac delta 
coefficient is so strong that the changes of regularity (i.e. appearance of the facets) happen instantly. 

We will state our results in the Section below, the proofs will be presented in the further Sections. 
Here, we present the outline of the rest of the paper. We show the existence of weak solutions in 
Section [2l uniqueness is the content of that Section, too. The qualitative analysis is based on the 
semi-discretization which is performed in Section [3] Our goal is to make some of the properties more 
apparent. Namely, we want to show that facets (i.e. intervals where tp = A s + s has a constant value 
equal to one of the aj's) form instantaneously. In Section [5j we show further geometric properties of 
solutions, namely the curve becomes convex (i.e. the angle ip becomes monotone) in finite time. In 
addition, we show that solutions become fully faceted in finite time, i.e. the solution is composed only 
of facets. These two events are not correlated in time. Finally, we show that our solutions converge to 
a special solution which we call minimal. 
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2 The main results 



Here, we present our results. We begin by noticing that if J is given by d 1 -4b . then the meaning of 
dl.21 ) is not clear at all because its right-hand-side formally becomes 

3 

A* = |^5 fc |_^(s + A s )A ss . 

k=0 

Hence the above equation can be viewed as a generalization of equation (11.5I ). 

We will use the tools of the convex analysis to interpret it. Due to convexity of J its subdifferential 
is always well-defined. Since in general d^J{4>) is not a singleton it is necessary to find its proper 
selection, in particular (11.21) takes the form, 

A* € J^V(A s + s), inSx(0,T), 

A(s,0) = A (s), on S, 

A(27r,i) = A(0,i), fort>0, 

where S is the unit circle. 

In other words, we have to find (weakly) differentiable selections of J(A S + s). Thus, we are 
lead to the following notion of a weak solution to (1 1 - 2b - 

Definition 2.1. We say that A G C([0,T\; L 2 (S)), such that A s G L^O, T; TV(S)) is a weak 
solution to f !2. 1 I ), if there exists a function $7 G Li(0, T; W{(S)) such that £l(s, t) G dI(A s + s) a.e., 
and for any function h in C°°(S) it holds 



f A t h = - f{n-s)h s + [ h. 
Js Js Js 



IS JS JS 

With this definition we can show the following existence result. 

Theorem 2.1. Let us suppose that J is defined by[L4\). A € Li(S) and A , s G TV(S) , then there 
exists A G C a (0, T; L 2 {S)) with a > 0, additionally 

A s G Lqo (0, T; TV (5") ) and A t G L 2 (0, T; L 2 (5)) 

such that it is a unique weak solution to j2.1\) . 



The proof will be achieved through an approximation procedure, it is performed in Section 3. 
Moreover, we show uniqueness of the solution constructed here, this is the content of Theorem l3.1l in 
Section 3. 

However, our main goal is to describe precisely qualitative properties of solutions to (I2.ll ). As a 
motivation, we present a special type of solutions, which we will call minimal solutions, which are 
given explicitly, one of them is given here, (see also §3.2), 



A(s,t) = I (p(u)du + t, 

where 



o 



7T 3tt Stt 7tt 

= 4*[o,|)( s ) + -jX\l^){s) + -jXfr*){s) + -jX { ^,2k)( s )- ( 2 - 2 ) 
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It is a matter of an easy exercise to see that A defined above with Q(x,t) = x is indeed a weak solution 
to (I2.ll ). In fact, this an asymptotic profile, which can be reached in finite time. 

We will keep in mind this example while developing the proper class of regular solution. The idea 
is that we want to extend properties of convex solutions to a more general class, hence we introduce a 
class of J-regular function, where restrictions on regularity depend on function J from (1 1 -4b - 

Firstly, we define the space of functions which are helpful to describe the regularity of the deriva- 
tive of our solutions. We recall that any function <fi G TV is a difference of two monotone functions. 
Thus, we shall call a multifunction (ft : [0, 2ir) — ► 2 R a maximal TV function if it is a difference of 
two maximal monotone multifunctions and one of them is continuous. 

Definition 2.2. We say that a maximal TV multivalued function <p : [0, 2ir) — ► M is J-regular, i.e. 
4> S J-R[0, 2tt), provided that the set 

H(0) = {s € [0, 2vr) : <£(s) 3 a k for k = 0,1, 2, 3} 

consists of a finite number of connected components, i.e. we allow only isolated intervals or iso- 
lated points. Additionally, on any connected subset [0, 2-7r) \ E function <j) takes its values in interval 
«fc + § ) for some k = 0, . . . , 3, modulo 2n - see (1.3\) . 

For each <ft 6 J-R[0, 2tt) we define a function K : J-R[0, 2ir) — > N by the formula 
K{4>) = the number of connected components of the set H(</>). 
Additionally we put 

||0||j-R[O,27r) = ll0l|TV[O,27r) +K (</)). 

Let us note that the J-R class does not form a Banach space. It is not a linear space. In order to 
formulate the meaning of solutions, first we define the composition of J-R functions with dJ. Because 
of the complex structure the definition is long. 

Definition 2.3. We define the composition dJoA, 

dJoA:[a,b] -» [e,f], 

where A : [a, b] — ► [c, d] is an J-R function and dJ : [c, d] — > [e, /] as follows: 

To begin with, we decompose the domain [a, b] into three disjoint parts [a, b] = V r UDjU T> s , 
where 

V s = {s <E [a, b] : A(s) = [c s ,d s ] and c s < d s }; 

(2.3) 

v f = {Ufc(°fcA) : A \(a k ,b k ) = c fc> where °k is a constant}; V r = [a,b] \ (V s UV f ). 

Then, the composition is defined in three steps: 

1. For each s 6 D r the set A(s) is a singleton, thus the composition is given in the classical way 

dJoA{s) = dJ{A{s)) for s£V r . (2.4) 

2. In the case s G T>t the definition is "unnatural". For a given set (aj;,6fc) C T>t we have 
A\r akj l, k \ = Cfe. If dJ(ck) is single-valued, then for s € (ak, bk) we have, 

dJoA(s) = &c k )}. 



6 



However, if dJ(c k ) is multivalued, i.e. dJ(ck) = \a k ,j3 k \, then the definition is not immediate. 
We have to consider four cases related to the behavior of multifunction A in a neighborhood of interval 
(a k , b k ). The regularity properties of the J-R class imply the necessity to consider the following four 
cases (for small e > 0): 

(i) A is increasing, i.e. A(s) < c k for s G (a k — e, a k ) and A(s) > c k for s G 6^ + e); 

(ii) A is decreasing, i.e. > for s G (a& — e, a^) and A(s) < for s G 6fc + e); 

(iii) A is convex, i.e. A(s) > c k for s G (a^ — e, a*) and > for s G (6fc, + e); 

(iv) A is concave, i.e. A(s) < c k for s G (a^ — e, a^) and < c k for s G (6&, b k + e). 
The case (i) we put 

dJoA(t) = x k (t - b k ) + y k (t - a k ) for t G (a k , b k ), (2.5) 

where = °\ and = , /3fc . 
For case (ii) we put 

dJoA(t) = x k (t -b k ) + y k (t - a k ) for i G (a k , b k ), (2.6) 

where x k = and % = , " fc . 

When we deal with case (iii) we set 

8JoA{t) = k for t G (a k , b k ). (2.7) 

Finally, if (iv) holds, then we put 

dJoA(t) = a k for t G (a k ,b k ). (2.8) 

3. In the last case, if s G V s our definition is just a consequence of first two steps. Since set 
V s consists of a countable number of points we consider each of them separately. We have A(d k ) = 
[e k ,fk] withe fc / f k , then 

SJoA(dfc) = [limsup dJoA(t),limmf d Jo A(t)]. (2.9) 

Definition 12.3 1 is complete. 

Thanks to the J-R regularity of A, the above limits are well defined. As a result, we are able to 
omit point from V s in (I2.3I ). We note that the above construction guarantees that 

dJoA : [a, b] — > [e, /] is a J-R function. 

After having completed the definition we make additional comments on step 2. Formulae (12.51 )- 
(12.81 ) are immediate consequences of the pointwise approximation of the considered function by 
smooth functions. The presented composition agrees with the results from [MRy ], where a stationary 
version of the problem has been considered. In particular, our definition follows from a requirement: 
if A is maximal monotone then we expect 

A~ x oA = Id. 

Moreover, the composition of two maximal increasing functions is maximal increasing. Another 
point, which should be emphasized, is the nonlocal character of the above definition. Step 3 depends 
on step 2, so steps 1 and 2 should be performed at the very beginning. 

Now we are prepared to introduce the main definition. 
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Definition 2.4. We say that a function A : S ■ ?i is an almost classical solution to system di.2D iff A 
is a weak solution with Q = dJo[A s + s], A s + s G Loo(0, T; J-R[0, 27r)) amc? 

A* = £dJo[A s + s] in [S x ((0,T) \ AT)] \ [j 0<t<T dE(A s (-,t) + s) x {t}, 
A\ t= o = A on S, 

where N is finite and dE denotes the boundary of set E. 

The main point of Definition 12.31 is to determine the composition appearing on the RHS of the 
equation on sets, where the solution and d J are singular. Note that equation ()2.10|) i is fulfilled in 
the classical sense except for finite number of point for each t G (0, T) \ N. This is so due to the 
definition of set 3(A s (-,i) + s) implying that its boundary consists of finite number of points. It is 
easy to see that the minimal solutions (12.21 ) fulfills Definition 12.41 

The main result of our considerations is the following. 

Theorem 2.2. Let Ao he such that Aq iS + s G J-R[0, 2ir), then there exists a unique almost classical 
solution to system (2.1b conforming to Definition \2.4[ 

In fact this is a statement about the regularity of weak solutions. Theorem 12.21 is a result of the 
semi-discretization of system (11.21) . At this level, we will be able to show that facets must appear, 
as suggested by the pictures in the Introduction. The semi-discretization will determine the RHS of 
(12.101 ) on sets where the solution falls into the singular part of d J. We will obtain that on these sets 
the term d J is constant on each connected part (or time dependent for the evolutionary system). Next, 
by the elementary means we will show that the semi-discretization tends uniformly to the solutions 
obtained by Theorem l2.ll However, performing a rigorous proof that we indeed constructed an almost 
normal solution requires more work on the structure of weak solutions, which is the content of Section 
[5J Thus, it will be postponed until the end of this part. 

At the end, in Section [5j we deeply go into the qualitative analysis of the evolution showing the 
convexincation effect and convergence to the minimal solutions. Since we know that facets must 
appear and the solutions are unique we are in a position to construct quite explicit solutions. We 
are able to follow their qualitative changes. This is made precise in Theorem 15.11 In particular we 
show instantaneous creation of facets. For the sake of this study we show a comparison principle in 
subsection [5] 1 . Moreover, we show that the evolution of facets is governed by a system of ODE's 
which are coupled if the facets interact, this is explained in Section [5J A conclusion from our analysis 
is existence of a sequence of instances at which our solution gets simplified before it gets the final 
form of the asymptotic profile, i.e. the minimal solution. 



3 Existence of solutions 

In this Section we show an existence and uniqueness of weak solutions of (1 1 - 2b - We use the tools of 
the convex analysis to interpret it. In particular, we shall make the gradient flow structure of (11.21 ) 
transparent. However, the existence is shown by the method of regularization. Some of the statements 
are easier to interpret if they are written in the language of the 'angle' <p = A s + s. Here, ip plays 
the role of the angle between the normal to the curve and the xi-axis. Thus, for convex closed curves 
ip must be increasing, but we shall not require that, instead we admit ip being a functions of bounded 
total variations, i.e., <p(',t) G TV(S), in particular ip G L°°(S) and it may be discontinuous though. 
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3.1 The proof of the general existence result 

We present a proof of our existence result, Theorem l2.ll It will be achieved through an approximation 
procedure. For any e > we set 



where p e is a standard mollifier kernel, with support in (— e, e). Let us note properties of the approxi- 
mation J £ : 

(a) J e e C°°(M); 

(b) J 6 is strictly monotone; 



(d) -^J e {x) - ex = -^J(x) for x such that \x - > e for k = 0, 1, 2, 3. 
We start with existence of the regularized system. 

Lemma 3.1. Let us suppose that J e is defined by O.il) and A is smooth and 2tt -periodic. Then, for 
any T > there exists a unique, smooth solution to the regularized problem, 

A e t = ^- s -^J e (K + s), inSx(0,T), 

A e (s,0) = A (s), onS, (3.2) 
A e (s + 2ir,t) = A e (s,t), for t > 0. 

Proof. By properties (a), (b) (c) and (d) of J e , see d3.ll ), the existence and uniqueness of smooth 
solutions to (I3.2I ). is guaranteed by the standard theory of parabolic systems, see BLSUL □ 
We now study properties of established solutions. 

Lemma 3.2. Let us suppose that A e is a smooth solution to 0.21) . 

(a) If for B,ieM and the initial datum satisfies a < (Aq s (s) + s) < b, then, for allt <T we have 




(3.1) 



a < (A e s (s,t) + s) < b. 



(b) If moreover, (Aq )S (s) + s) s G Li(0, 27r), then, for all t <T we have 



(A e s (s,t) + a), e Loo (0,^^(0,2^)). 



Proof. We use the maximum principle. First of all, we differentiate (13.21 ) with respect to s, 




We notice (s + A e s ) t = A € st . We set w 



(s + A e s ), hence we obtain the equation for w, 



w t = —(a{s,t)w s ) 



(3.3) 
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where by (13.11 ) we have a(s,t) = (s + A e s ) > e > 0. Hence, by the maximum principle we obtain 
(a). 

To prove (b) we note that from (13.31 ) we obtain 

d 2 

w s t = j^(a(s,t)w s ). (3.4) 

By Lemma 13.11 our solutions are smooth. In oder to finish the proof of (b) it is enough to integrate 
(13.41 ) over sets {w s > 0} and {w s < 0} to reach, 

w s dx < and — / w s dx > 0. (3.5) 



dt J{ Ws> o} dt J{ Ws< q} 

□ 

Having established this Lemma, we will obtain estimates for the spatial derivative of solution 

A. 

Corollary 3.1. There is a constant M independent of e and T such that 

II^ € |ILoo(5x(0,T)) < M, ll ( / ; ' e (-)*)llL oc (0,T;Ty[0,27r)) < M - 

Proof. The first part follows from Lemma l3T2l (a) directly, because ip e = A e s + s. The second part is the 
result of Lemma l3T2l (b). combined with the properties of approximation of TV functions in L\. □ 
We want to show that the estimates for A e will persist after passing to the limit with e. 

Lemma 3.3. Let us suppose that A e converges weakly in L 2 (S x (0, T)) to A. If (A e s + s) s > in 
V'{S), then (A s + s) s > as well in V'(S). 

Proof. Indeed, if h € T>(S) is positive, then < J s {A e s + s)h s . The inequality holds after taking the 
Umit. □ 

Lemma 3.4. There is a constant independent of e such that 

\ (A e ) 2 dxdt <M, / / [{A%) 2 + (A t e ) 2 ] dxdt < M. 
ii Jo Jo Jo 

Proof. The bound on Jq (A 6 ) 2 is trivial, due to L°° estimates established in previous lemmas. 

,T ,2 

h Jo 



T i-2-k i-T r 2-n 

(A e t ) 2 dsdt 

JO 



Similarly, the bounds in Corollary 13. II imply that J* Q T ^ {A e x ) 2 < M. We shall calculate the last 
integral with the help of integration by parts, 

£ £ n Kt-^nK + s) dsdt + £ K^J € (K + sWsZl 71 dt 

T J e (M*))ds- [ J € (<p(s,T))ds 
Js 

+ £ Al(0,t) U-r(At(0,t) + 2vr) - ±J*(Al(0,t))\ dt, 

Here, we also exploited periodicity of A. We notice that the difference ^J e (A|(0, t) + 2ir) 
^J e (A|(0,t)) equals exactly 2tt. Hence, 



T r-2-K r 2Tr 

(A e t ) 2 dsdt < / J e (M s )) ds + 27r(A e (0, T) - A e (0, 0)) < M 
o Jo Jo 
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due to the Corollary 13.11 □ 

Remark. We want to stress that the above estimate on A t is one of the most important differences 
between (11.11 ) and (|2~TT ). 

Now, we have enough information to select a weakly convergent subsequence, with properties 
announced in the theorem. 

Proposition 3.1. There exists a subsequence {e^} converging to zero, such that 

(a) A e * ->« A in W^(S x (0, T)); ipf -± p> s as measures in S x (0, T). 

(b) AeC([0,T),L 2 (S)). 

Proof. The first part of (a) is implied by Lemma l3~4l The second part of (a) follows from p e = A e s + s, 
and Lemmas [3 .21 13 .41 Part (b) follows from Lemma l3~72l and l3~74l and the embedding theorem (we have 
already proved A e G L 2 (0, T; W£(0, 2tt)) n W£(0, T; L x (0, 2tt))). □ 
The next step is to show that that the limit is indeed a solution. In particular, we have to pass 
to the limit in the non-linear term. First of all, we shall change the notation in order to make more 
transparent what we are doing. We want to find w(s,t) such that w s (s,t) = ip(s,i). By a simple 
integration of this formula and the definition of tp, we can see 

w(s,t) = ^s 2 + A(s,t), 

where we set w(0, t) = A(0,i). Hence, w s = <p and we can re-write the evolution problem as a 
gradient system 

w t £-^dJ(w s ), mSx(0,T), 
as 

w(s,0) = ^s 2 + A (s), for seS, (3.6) 
w(s,t) — is 2 is periodic for t e (0,T). 

If </?(•, 0) is increasing, then due to Lemma [3721 (b) and Lemma [3731 ip{-, t) is increasing as well, 
hence w(-,t) is convex. Obvious changes are required to write the system for the regularization 

w e (s,t) = ±s 2 + A e (s,t). 

Proposition 3.2. For any fixed t > and a sequence {e^} converging to zero there exists its subse- 
quence {ek} (not relabeled), such that for each x G [0, 2tt) the limit 

exists. Moreover, Q(x,t) G dJ((p(x,t)) for almost every x € [0,27r). 

Remark. It is important for us to make the selection of the subsequence independently of t. 

Proof. Indeed, once we fix t > 0, we may recall that p e (-,t) G TV as well as ■£^J € (tp t )(x,t) G 
TV. Hence, by Helly's convergence theorem there exists a subsequence such that these sequences 
converge. Using the new notation, we write, 

lim p e (x, t) = w x (x,t), lim — J e (w e (x,t)) = Q(x,t). 
e^O e^o dp 
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Now, we shall show that for each point x the number U(x, t) belongs to dJ(w x (x, t)). Since the 
functions J e are convex, we have the inequality 



2-7T r2ir J 

J e (w x (x, t) + h x {x)) — J e (w x (x, t)) dx > / — J € (w x (x,t))h x (x)dx, 
o Jo 



dip 



for each h G Cg°(0, 2ir). We know that w e and -^J e {w%(x,t)) have pointwise limits, which are 
bounded, hence after passing to limit our claim will follow, 

2tt r-2-K 

J(w x (x,t) + h x (x)) — J(w x (x,t)) dx > / Q(x, t)h x (x) dx. □ 
o Jo 

We finish the proof of Theorem \2.1\ By previous Lemmas there exists a sequence A e which 
converges weakly in W%(S x (0, T)). In particular, if h G C£°(0, 2vr), t > and r > is arbitrary, 
then we see 

/ + / Afhdsdt' = f + [ ^--^-JHAl + s)/idscft' = - f + [ ^- J e (K\ + s)h s dsdt'. 
Jt-r Js Jt-r Js ds dip J t _ T J s dip 

Since -g- J e (A € s + s) is bounded, it converges weak-* in ^((O, 2ir) x (0, T)) to 0. We have to show 
that Q(s, t) G dJ(A s + s). First we notice that we may pass to the limit in the above integral identity, 

t+r i- rt+r i- 

/ A t (s,t')h(s)dsdt' = - / / n(s,t')h s (s) dsdt' . 

t-T JS Jt-T JS 

By the Lebesgue differentiation theorem we deduce, 

A t (s,t)h(s)ds = - / n(s,t)h 8 (s) ds (3.7) 
s Js 

for a.e. t G [0, T] for /i G W^S 1 ) (we used the fact that is not distinguished on S). In principle, the 
set G = {t G [0,T] : (|377T > holds} depends upon /i, i.e. G = We shall see, that in fact we 

can choose G independently of h. Let us recall that W^iS) is separable and let us suppose that D is 
a dense, countable subset of W^iS). Of course, Q = HheD i s a set °f f un measure. Let us then 
take t G £7 and /i G C°°(S'). Let us suppose that {h n } is a sequence in C°°(S') converging to h in the 
W2(£)-nann. Then, 



A t (s,t)h n (s)ds = - / Q(s,t)(h n ) s (s)ds 
s Js 

for all teG. We may pass to the limit with n on both sides, thus we reach, 

/ A t (s,t)h(s)ds = - / fl(s,t)h s (s)ds. 
Js Js 

In other words, (ED holds for all /i G C°°(S') and all i G Q. 

If we now fix t G ^, we next apply Proposition 13.21 to deduce that 0(s,i) G dJ(A s (s,t) + s). 
Hence the limit, A, is indeed a weak solution. □ 

Now, we are going to prove uniqueness. 

Theorem 3.1. If A\ i = 1,2 are two solutions with A 1 (s,0) = A 2 (s,0), then A^s,*) = A 2 (s,t), 
fort < T. 
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Proof. If A*, % = 1, 2, are weak solutions, then by the definition of weak solutions we have 

/ h\hds = - / (fi* - s)h s ds+ / hds, 
Js Js Js 

where w % G — dJ and h is in i? 1 . We subtract these two identities for A 2 and A 1 , then we take 
(A 1 — A 2 ) as a the test function. Finally, the integration over (0, t),t < T yields 

f \ Jt-( a1 " a2 ) 2 dsdt = ~ I \ ( Ql ~ ^ 2 )( Al " a2 )^ dsdL 
Jo Js 2 dt J J s 

Monotonicity of d J implies that ^A 1 - A 2 || 2 2(s ^(t) < 0. Hence, ((A 1 - A 2 || 2 2 ^(i) = for any 
t < T. □ 

3.2 Minimal solutions 

It is well-known that important information about the studied system is provided by special solutions, 
like traveling waves, self-similar solutions and other symmetry solutions. We can not talk about self- 
similar solutions because our systems lacks direct geometrical interpretation, however we may look 
for special ones, which we named minimal solutions. 

In the theory of curvature flows it is natural to anticipate existence of curves such that their cur- 
vature is constant, but may change in time. Here, we ask if there exists such a solution (p to (11.21 ) 
that 

-rdJUp) B k, hence dJUp) 3 ks + s*, (3.8) 

as 

where s* is appropriately chosen, e.g. s* = ~ and \k\ = 1. The last restriction is of geometric nature, 
namely we want that for any a G R the image of S by d J be contained in an interval no longer than 
2vr. 

In fact, we may come up with explicit formulas. One for k = 1 is provided by formula (12.21) . It is 
then obvious that 

<p{s):={dJ)- l {s + s*), (3.9) 

as in |MRy| and in Section 2. Moreover, <p(s) ds = 2tt 2 = sds. By the reversal of the 
orientation, we immediately obtain the solution for k = — 1, 

yyr 57T 371" IT 

We choose <p(s) := (p\ (s), which is given by (12.21) . because we prefer to have (p an increasing function. 

As a result, A defined by A(s, t) = J* Q S <p(u) du+F(t) is indeed 2ir periodic in s and it is a solution 
to (1 1 -2b - Here, we must take F(t) = A + t. One can check in a straightforward manner that indeed A 
solves (|3.9b . This is indeed so, because we have found tp(s) = A s (s) + s and $7 is a section of dl{(p), 
namely, Q(s,t) = s, which satisfies (I3.9I ). If we take A = 0, then A satisfies the initial condition: 
A(s,0) = f s <p(u)du. 

4 The semi-discretization 

In this part we examine the semi-discretization of (I2.ll ). Our goals are not only to establish existence 
for the presented scheme, but also to show qualitative properties of the obtained solutions. In partic- 
ular our considerations will explain the appearance of facets. Finally, we prove the convergence of 
solutions of the semi-discretization to the solutions obtained in Section 3. 
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We define the semi-discretization in time of system (12.11) as follows 

A£(8) "^" 1(a) €|aj[AtW + a ] (4.D 

and A£(0) = A^(27r) and (A°) s = (f> for = 1, . . . , [T//i]; or equivalently equation (|4~TT) can be 
stated 

A fe( s )_^j[ A fe js(s) + s]9A fc-l( s ). (4.2) 
We establish existence of solution to this problem. 

Lemma 4.1. Let us suppose that an absolutely continuous function v is such that v s = (p G 

ry[0, 27r), then there exists u G AC([0, 27r)) such that u s G TV", which is a solutions to ||4JD . 
i.e. 

u — v £ h^dJ(u s ) (4.3) 
as 

with u(0) + ^(27r) 2 = u(27r) and the following bound is valid 

\\u s \\tv < \\v s \\tv- (4-4) 



Remark. Our understanding of (14.31 ) is the same as that of (12.11 ). i.e., there exists uj G Wj ([0, 27r)), 
such that u;(x) G dJ(u s ) and it — u = h-^uj. 

We also note that u and u appearing in this Lemma need not be periodic, on the other hand A(-, t) 
and A^(-) are periodic. 

Proof. Let us notice that if u is a solution to ( 14.31 ), then belongs to the subdifferential of the functional 

J(u)= [hJ{u s ) + -(u-v) 2 }, for u G AC([0, 2ir), u s G TV. 
Jo 2 

i.e. u is a minimizer of J. To be precise, we define J on L 2 (S) by the above formula for u G 
^4C([0, 2))7r with u s G and we put J{u) = +oo for u belonging to the complement of this set. 
In order to solve ( 14.31) , we consider a family of regularized problems, 

J e (u)= [ [hJ e (u s ) + l(u-v) 2 }, 

Jo 1 

where J e is the same regularization of J that we used in d3.ll ). 

The functional J e is well-defined, convex and coercive on the standard Sobolev space Wi(0, 2tt), 
thus it possesses a unique minimizer u e . Now, we apply again the methods used in Section 3.1 to 
show existence of a weak solution of the evolution problem (1 1 -2b - The regularization of system (14.31) 
leads to the following equation 

d 2 8T 

Uss ~ Is 2 mp 2 s ~ Vss ' 

By repeating the argument for (13.41 ). we get < H^Lll^i- Passing to the limit with e — > 

yields (t4~4l >. 

In addition we have the following bounds f 2lT (u e ) 2 dx < M, f 27T (u e x ) 2 dx < M. In order to 
prove them we follow the lines of reasoning of Corollary 13. II and Lemma |3~4l These bounds suffice 
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to show existence of a subsequence {e^} converging to zero, such that 
(a) u ek — 1 u in W% (0, 2-7r); u e s k s — 1 u ss as measures. 

Subsequently, by Helly's theorem we conclude existence of the pointwise limits (for another sub- 
sequence {efc}, not relabeled) 

lim tp e (x) = u x (x), lim — J e (u|(x)) = O(x). 

t—>oo t^co dip 

Moreover, Q(x) G dJ(u x (x)) for each x G [0, 2tt). 

Now, we show uniqueness of solutions, constructed in Lemma |4~T1 

Lemma 4.2. Let v G AC([0, 2tt)), t> s G TV([0, 27r)), tnen tnere exists at most one weak solution 
u G AC([0,27r)), m s G rV(5) to probJem ifO) 



Proof. Let us suppose that there are two solutions to (14.31 ). u\ i = 1,2. By the definition, there 
are two functions u>j G dJ(u l s ), i = 1,2, such that 

u — v = h—oJi, % = 1, 2. 
as 

After subtracting these two equations and multiplying them by m — U2 and integrating over [0, 2tt) 

f / ^ ^ \/1 9\ 

/ cji — cJ2)(n — it ) ds = 0. 

Jo «s ds 



we see 

f27T 

I 1 2 M 2 
\U — U \\ 



The integration by parts leads us to 

r-2n 

= Wu 1 - u 2 \\ 2 + \ h(ui - uj 2 ){ul - u 2 ) ds > Wu 1 - u 2 \\ 2 > 0. 
Jo 

As a result u 1 = u 2 . □ 
In order to finish our preparations, we introduce the sets of preferred orientation which dominate 
the behavior of solutions. Let us suppose, that w is absolutely continuous and w s G TV, then at any 
point s, the left derivative w~, as well as the right derivative wf are well-defined, hence we may set 

dw(s) = {twJ + (1 - r)wt : r G [0, 1]}. (4.5) 

If w is convex, then dw is the well-known subdifferential of w. 
Now, for each I = 0, 1, 2, 3, we set 

Ei(w s ) = {s G [0, 27r] : it) is differentiable at s and u> s (s) = or G 5w(s)} (4.6) 

Furthermore, we set S(w a ) = Uf=o —iC" 7 *)- 

The result, delivering the main properties of solutions, is the following. 

Theorem 4.1. Let 4>q = A° + s G J-R[0, 7r). Then a solution {A^} to problem if4.il) exists, it is 
unique and it satisfies the following bound 

Pfc,s + S l|j-R[0,27r) < ||A ,s + s||j-R[0,27r)- ( 4 - 7 ) 

Moreover, we have 

H(Ati )S + a) C H(At + s) and K(A* s + s) < K^ 1 + s) (4.8) 
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and 

sup sup |Sf \ Sf _1 | < C(Vf>) + /i 1/2 ). (4.9) 

fe 2=0,1,2,3 

where V(s) — > as s — > and V is determined by the initial datum <pQ. Moreover, on connected 
components of the set E^ 1 \ ((jf=o \ S f _1 ) 

^-5J[A^ s + s] is constant. (4.10) 

Proof. By Lemmas 14.11 and 14.21 we conclude existence of the sequence of solutions to the semi- 
discretization, the solutions are such that A^ belong to TV(S). It is enough to restate the equation 
(14.21 ) as follows: 

u-h^-dJ\u s ]=v (4.11) 
ds 

with u = X ! l + \s 2 and v = A^ _1 + \s 2 , and boundary condition u(0) + 2ir 2 = u(2n). 

The set, where function J[u s ] is singular, i.e. 3(it a ), plays the key role. Our first task is to prove 
the inclusion from (14.8b - Note that in a neighborhood of any point s ^ E(u s ) function 9J[u s (-)] is 
constant, hence we get u(s) = v(s). Thus, we point the first feature of solutions to (14.111) 

u(s) = v(s) for s £ (0, 2vr) \ E(u s ). (4.12) 

From (14.121) we deduce that if s £ E(v s ), then s E(u s ). Subsequently, we get E(y s ) C E(u s ) 
which proves the inclusion from (14- 8b - Thus, the isolated elements stay isolated or merge with other 
elements. From this we obtain that K(v s ) > K{u s ) what ends the proof of line (14.8b - 

By properties (14.121) . (14.81 ) and the estimate from Lemma (14.11) . we immediately deduce estimate 
(14.7b - In particular, what we gain is a uniform bound in L OQ (S) on {A^ s }. 

The set E(u s ) is defined as the sum of Uf =0 E[(u s ), thus without loss of generality we can con- 
centrate our attention on one of them, e.g. on the set ^(u^) - see (14.6b . From the J-R-regularity of 
u s set S2 (it s ) is a sum of closed intervals, so we take one of them, say, 

7T 

[a_,a+] C E 2 (it s ) and n s | (a _ ja+) = -. (4.13) 

Recalling the required regularity of the functions in the J-R-class, we find e > such that one of the 
four following possibilities holds: 

W ^s|(a-— e,a_) 4; u s I (a+,a++e) < 4~j 

(*'*) u s|(a_-e,a_) < f> u s\(a + ,a++e) > f ! ^ 

(m) lts|(a_— e,a_) ^ 4; u s I (a+,a++e) > 4") 

O) «s|(a_-e,a-) < f> ^l(a+,a++e) < f • 

Subsequently, we integrate (14.11b over (a_ — e, a + + e) to get 

/ uds - h{dJ[u s ]\l + _t\) = / uds. (4.15) 

J CL——e Ja——e 

After passing with e — > + , we obtain - according to the above four cases ( 14. lib - the following 
identities 

(i) uds — /i| = vds (convexity), 

(ii) uds + hr; = fa vds (concavity), (4.16) 
(Hi) and (iv) f° itcZs = J" vds (monotonicity). 
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In our present analysis, we essentially use the fact that the energy density function J is defined 
by a square. Due to the definition of J, see (11.41 ). formula (14.161 ) exhausts all the possibilities of the 
behavior of u s . For more complex polygons, we would have to discuss more possible types of facets 
- here, there are just four of them. 

We keep considering the interval [a_, a+] C E 2 (u s ), see (14. 13b - Let us introduce a set 

n = ([a_,a + ]nH( I)s ))\(H l) ( ? i s )UH 1 (t ls )UH3( Us )) 1 (4.17) 
then by the properties of sets H, we deduce that 

(u — v)\u = Ch is constant. (4.18) 
The sign of constant Ch is determined by the geometrical properties of cases in (14- 16b - We have 

C h > for (i), C h < for (it) and C h = for (in) and (iv). (4.19) 
Also identity (14.181 ) and equation (14.111 ) yield 



^-dJ[u s 
as 



Ch a d a t\ 
= — — and —oJ\u s 

n h ds 



= 0. (4.20) 

(0,2tt)\H(m s ) 



Thus, we proved ( I4.10I ). 

Next, we are going to study ( 14.91 ). From the analysis of (14.1 II) . we conclude that 

\W-v\\ Ll{s) <h^K{4> ). (4.21) 

Additionally, from (14.81 ) we have also that u, v 6 W^ Q (S), thus simple considerations lead us to the 
following bound 

\W-v\\ Loc{s) </i 1/2 C(0o). (4.22) 
In order to measure the set Ei(u s ) \ Ej(t> s ) we split it into two parts 

E^us) \ Et(v s ) = [(S,(« a ) \ Ei(v a )) n S(t;.)] U (H,(«.) \ E(v s )) = Hi U n 2 . (4.23) 

Let us consider IIi. On this set we watch the evolution of the intersection of facets. Thanks to the 
full information about the direction of this facet, we deduce immediately that 

III! | < C{4>v)h 1 ' 2 (4.24) 

The number of possible intersections is controlled by K(4>o). 

To estimate 1I2, let us note that this set is a subset of E(Xq :S + s), thus in the general case we can 
say only 

|II 2 | < V(h), (4.25) 

where V(s) — > as s — > and V is determined by the initial datum. Assuming strict convexity of 
initial domain we would obtain V(h) ~ /1 1 / 3 - see the example at the end of subsection [5]2. □ 
Theorem 14. H is proved. 

Next, we show that sequences {A^} converge to solutions of the original problem. We will com- 
pare solutions given by Theorem 12.11 and Lemma 14.11 in particular, all assumptions of Theorem 12.11 
are not required. We follow the standard procedure which is valid for parabolic operators (see IMRall). 

Our next task is to show the following lemma. 
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Lemma 4.3. Let A and {X^,} be solutions to problems (f2.il) and if4.il) respectively, then 



[T/h] 



\\Hs,t) - ^(a)X[fc,*+i)(*)IUi(o,r ; i a (S)) as h H 



(4.26) 



fc=0 



If the initiai datum fulfills the assumptions of Theorem \2.1[ i.e. Aq )S € TV(5), then 



[T/h] 



\\A(s,t)- ^h(s)X[k,k + i)(t)\\ 



k=0 



L p (0,T;W?- e (S)) 



as h -> H 



(4.27) 



for any 1 < p < oo and e > 0. 



Proof. From the properties of solutions to problem (I2.ll ). we know that At £ £2(0, T; L,2(S)). It 
follows that 



A(s,i) - A(s,i- h) 



A t («,t) 



as /j -> 0^ 



(4.28) 



Li(h,T;L 2 (S)) 



For fixed /i > we denote 



Rh(s,t) = A t (s,i), 



then the equation (I2.1l) i can be restated as follows 

A(s,i) - A(s,t- h) 



h 



-Trds 



0(s, i)7r s 4- Rh{s, t)irds 



(4.29) 



(4.30) 



for each 7r in C°°(S x (0, T)) and each selection f2(s, t) of multivalued function 5J[A s (s, i) + s] 
We want to compare the above system with the semi-discretization given in Section 0] 



h 



■irds 



u(s, t)irds 



where t € [kh, (k + l)h) and u>(s, t) is any section of dJ[X k s (s) + s]. 
Let us define 

A k (s,t) = A(s,t)-X k h (s,t), 
provided t E [kh, (k + l)h), then from ( 14.301 ) and ( 14.3 lb we deduce 



{0(s, f.) — lo(s, t)} tt s + Rh(s, thirds. 



is h 
Taking in (14.331) as a test function A k (t, s), we get 

W) \A k (s,t)\ 2 ds = J (0 27r) A^s.tJA*" 1 ^,* - fc)cfe 
-^W) («(a,t) - w(s,t)) (A a (a,t) - A*,(M)) 
+/l /(o,2^) Rh{s,t)A k (s,t)ds, 



(4.31) 



(4.32) 



(4.33) 



(4.34) 
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but the monotonicity of d J implies 



f t) - u(s, t)) (a s {s, t) - \ k s (s, t) 



ds > 0. (4.35) 



So, defining a k (t) = \\A k (-,t)\\ L2 ( S ), by the Schwarz inequality, we get from ( 14.291 ) the following 
inequality 

a k (t) < a k -\t- h) + hrl(t) for t G [kh, (k + l)h), (4.36) 
where r k h (t) = | \R h (; t)\ \ L2i0 , 2n y Thus <E36]> yields 

k 

a k (t) < a°(t- kh) + ^2hr l h (t- (k - l)h) with t G [kh, (k + l)h). (4.37) 
i=i 

Integrating (14.371 ) over i£ [kh, {k + l)h), we get 

r(k+l)h j-h f-T 

/ a k (t)dt< oP(r)dT + h \\R h (;t)\\ L2{0>2w) dt (4.38) 

Jkh JO JO 

for T > (k + l)h. Introducing function a(t) = J2i=o a '(*)X[«/i,(«+i)/i)( <: ) with L = [T/h], from 
( 14.381 ) we get 

/ a(t)dt < /l 1/2 r||A t || L2(0iT;L2 ( ,27r)) +T\\R h \\ Ll (p,T;L 2 (0,2n)), (4-39) 
JO 

because the first term of the right-hand-side (RHS) of (14.391 ) is a consequence of the following estimate 

1 f h 1 f h f 4 

a°(t)dt<- \\ A t (-,T)dT\\ L2{QM) dt 



h Jn h 

,-h 



JO 



(4-40) 

Hence, from ( 14.391 ) and (14.281 ) we conclude ||ah||Li(o,T;L 2 (o,27r)) ^0 as h —>■ + and we get 
(14.26b - From the interpolation estimates and the results of Theorem l2.ll we deduce that for any p < oo 
and e > the convergence (14.271 ) is valid. Lemma 1431 is proved. □ 



5 Analysis of solutions 

The semi-discretization process proves that the set grows with k. One can show that the 

set |Ji>o ^( u h,s{', t)) may be estimated from below to show that it survives the limiting process as 
h — ► 0. This may be achieved by the analysis of the semi-discretization procedure, but this seems 
tedious. We propose an alternative approach by the construction of an explicit solution to dl.2| ) for 
data in (po G J-R. By uniqueness result, see Theorem 12. 1[ this is the solution. 

We shall assume in this Section that <p = w s belongs to J-R and this is the case for the initial data 
ifo of system (12.1b . As a result of the definition of the J-R class we see 

3 N 
1=0 k=l 
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where 

C<£ and £<£^i> k = l,...,N , (5.2) 

(with the understanding £/v +i = ?F + ^vr). Moreover, each interval [£;7,£^~] is a connected com- 
ponent of one of the sets Ei(ipo), I = 0, 1, 2, 3. We shall also adopt the convention that < £r and 
possibly > 2vr, but ^ - 2vr < £7. 

If [£7,£7] is one of the connected components of Hj(y), then we will call by a. facet the set 

^ = Fk(£k,tk) = i( x >y) e R2 : y = a e Kfc.^ib"]}- The interval [£7,£+] will be called 

the pre-image of facet F^. Let us stress that we admit £7 = i.e. a facet degenerated to a point as 
well as £7, = £7, i.e. we expect interaction of facets. We shall see that the generic initial data lead 
to the facet creation (from the degenerate ones) and their interaction. We show that facets are formed 
instantaneously from the data. At this point we mention that creation of interacting facets leads to 
additional difficulties and this process is handled separately. 

We will come up with an explicit formula. Once we check that indeed this formula yields a 
solution to equation (12.11) . we will be assured that this is the unique solution we seek. Subsequently, 
we shall see that solutions get convexified, i.e. after some finite time the angle becomes increasing, 
hence w becomes convex. Finally, we study interaction of facets. We will prove that w(-,t) becomes 
a minimal solution at the limit time. 

It will be also convenient to say that a facet Ff.(£7 , £iT), has zero curvature, if [£~,£iT] is a 
connected component of H(</>) and there exists an open interval (A,B), containing [£7,£jjT] such 
that w s is not monotone on any interval (a, b), satisfying 

[^4 + ]c(a,fe)c(AB). 

Furthermore, we say that a facet F^ = i*fc(£7 , fit) is regular if £7 < £+. Otherwise, we say that 
Fk is degenerate. If w s € J-R is such that the graph of w contains degenerate facets, then we say that 
facets are created in solutions to (12.11) . 

Finally, we say that facets Fi, . . . , Fi +r for r > 0, interact (or are interacting) if Fk n Fk+i, 
k = I, . . . ,1 + r — 1, is a singleton. We call a single facet F^ non-interacting, if it is not true that it 
interacts with any other facet. 

Thus, we have the total of eight combinations, we will treat each case separately. 

5.1 A comparison principle 

We are going to establish that solutions to equation (12.11 ) enjoy the expected comparison principle. 
This result is interesting for its own sake but also it is a useful tool analysis. We will apply it to show 
creation of interacting facets. 

We first recall the basic result (see, e.g. [S ]). 

Proposition 5.1. Let us suppose that u\, 112 are smooth solutions to a strongly parabolic equation 

u t = (a(x, u x )) x in S x (0, T) 

andu2(x,0) > ui(x,0), then U2(x,t) > u\(x,t) for all t G (0, T). □ 

With this result we may deduce the following comparison principle. 

Proposition 5.2. Let us suppose that K\, A2 are weak solutions to (12. J I) and Ai(x, 0) < &2(x, 0), 
then A 2 (x,t)> Ai (x, t) for all t £ (0, T) . 
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Proof. Since Ai(x, 0) < A 2 (x,0), we deduce that Af(x, 0) < A|(x,0), where k\, i = 1,2 are 
solutions to the regularized system (I3.2I ). Application of the preceding result yields 

A\{x,t) < A 2 (x,t). 

Since the point-wise limit exists we conclude that our proposition holds. □ 
We stress that no information about fij, i = 1, 2 is needed in the proof of the above result. 

5.2 Facet formation 

We shall see below that the evolution of a facet separated from other facets is governed by an ODE 
for its end-points, see (15.141 ) below. In the case of interacting facets their evolution is described by a 
system of ODE's d5J9l 

As we mentioned we admit facets F k degenerated to a single point at the initial instance to = 0. 
In this case the single ODE (15.1 lb and system ODE (15.191 ) become singular. While we can resolve 
satisfactorily the singularity of the single ODE, the analysis of the system is more difficult. In fact, 
we circumvent this problem by using the comparison principle to show creation of interacting facets. 

We shall use the notions and notation introduced above. In addition, in order to facilitate our 
construction we shall write 

x i-> a k (x - Sfc) + r k =: l k (x,s k ,T k ), 
where a k G A, s k G H (<p), T k G R. 

Theorem 5.1. Let us assume that <pq = wo.s G J-R and w is the unique solution to &2.1\) . We also 
assume that the set E(wo >s ) = Ufc=i[£fco> £fcol conditions i5. JD and (f5.2D . Then, there exists a 

finite sequence of time instances < to < H < • • • < < 00 an d & finite sequence of continuous 
functions 

£fc : Mi+i] -»R, i = 0,...t M -i, k = l,...Ni, 
^:[t M ,oo)^R, k=l,...N M = 4, 

where N > N x > . . . > N M = 4. 

The functions £!"(•)> £jT(") satisfying H5.2\) have the following properties: 

(a) £(0) = £ ; 

(b) < er(t) < < ^2 (*)<••• < ^ ^ 6" + 2vr, t G [U,U+i)i 

Ni 

(c) for t G [U,U+l) we have E((p(-,t)) = (J [^(t), £ (t)}, and each interval is a 

fc=l 

connected component of one of the sets H/(<^(-, t)), I = 0, 1, 2, 3. 

There exist functions r k : [tj, tj+i) — ► R, i = 0, . . . tu, k = 1, . . . iVj, and t^f+i = °°- They 
are such that the unique solution to (2.B with initial data ip(x, 0) = ipo(x) is given by the following 
formula for t G [£», tt+i)> i = 0, . . . , M 

Wx ^ = { «*>(*) ifx G [0,2tt) XUEaKfcW.tfWl (5 3) 

l ' J I lk(x,^(ti),T k (t))+w(^(U),ti) ifx €[£"(*),£+(*)], fc = l,...,^ 

Moreover, w x (-,t) is well-defined a.e., dw defined by (4.5\) belongs to J-R and 
[|£M-,i)||j-R< ||9™0 || J-R. 
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In addition, at each time instant ti, i = 0, . . . , M, one of the following happens: 

(i) One or more zero-curvature facets disappear, i.e. if one facet disappears atti, then 

< 4~ W < < ^ Q+1 (t), fort, < t < t l+l 

and 

I™ ^ -i(i) = lim ^ Q+1 (t) = 

where JS a subset of a connected component ofEi(ip(ti + i)), as a result Ni + \ < 

Ni. " 

(ii) One pair or more pairs of facets begin to interact, i.e. < (*) forti < t < ti+i and 

lim £_ 1 (t)=$_ 1 (t i+1 )=Z-(t i+1 )= lim £"(t). 
t— *t i+1 t-+t i+1 

The proof is achieved in a number of steps. Its major parts are separated as Lemmas. We start 
with constructing the £^'s. We first consider non-interaction during creation of facets, i.e. 

if£ = #, then Ci<C" n dC<C+r ( 5 - 4 ) 
However, the lemma below is valid without this restriction. 

Lemma 5.1. Let us suppose that w s = 92 G J-R and [£Z\£i!~] is a connected component ofBi(ip) 

and s k is its member. We assume that F k (^ , is not a zero curvature facet. 

(a) If g£ < £fc+i> then for sufficiently small r k of a proper sign, there exist ^ (r k ) such that 

=k($(T k ),$,T k ) + w($) and^(0) = £. (5.5) 



Moreover, the functions r k ^(ffc) are Lipschitz continuous, provided thatw s (t£) / ojfc. Other- 
wise, ^( r fe) 31-6 locally Lipschitz continuous. In addition, 

^$(T k ) = —-^ , = — 7 7^ fora.e.|r fc |G[0,e). (5.6) 

(b)If£+ < < = ^ ^+2 = - < £j+ r < ^+ r +P (in particular we admit 

er = & - 27r ^ ti3en 

w(C-i) +^-i(Cfc_ 1 (Tfc-i,r fe )),^+_ 1 ,r fe _i) = w((, + ) +l k (^(T k ,T k+1 )),£+,T k ) (5.7) 

for A: = 1 + 1, . . . , Z + r. Moreover, the functions {T k ,T k+ \) 1— ► £, k {r k , T k+ \), k = I + 1, . . . ,1 + r — 1 
are Lipschitz continuous. 

Proof. Before proceeding to the formal proof we will explain the situation by drawing a picture (where 
the subscript k is suppressed). The graph of w(-) and the line containing F(^ , £+) moved vertically 
by r intersect at x = £~ (r) and at x = £ + (t). 



(a) Since F k is not of zero curvature then by the fact that 9? € J-R it follows that w in a neigh- 
borhood of is either convex or concave. Let us consider the case of w being convex on 
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(a, b) D [£7 , the other case is similar. By convexity, any chord is above the graph of w. Thus, 
the line l k {-,£, k ,r k ) + f° r sufficiently small r k > intersects the graph of w at exactly two 

points, i.e. for r k > equation (15.51 ) has exactly two solutions. One of them, which is greater than £7 
is called £fc~( r fc)> the other one, smaller than £r is dubbed £7(7"fe)- The function 

x ^ w(x) - lk(x,Z+,Tk)-w($) =■ K( x ) ( 5 - 8 ) 
is increasing for x G [£7 , 6) and this interval is maximal with this property, while the function 

x i ^ io(x) - h(x,^,T k ) - w(£7) =: F7"(x) ( 5 - 9 ) 

and decreasing for x G (o, £7] and again this interval is maximal with this property. One can see this 
by taking the derivative of (15.81 ) and ( 15.91 ), because we have 



-j^( w ( x ) - h(x,^,T k )) = w'(x) -a k > w'(^) -a k >0 for a.e. x G 

and 

-fr,( w ( x ) ~ h(x,C^,Tk)) = w'(x) - a k < u/(£7) - a k < for a.e. x G (a, £7]. 

Thus, the function 6) G x i— > F^(x) (resp. (a, £7] i— > F7"(x) has a continuous inverse. As a 
result, for any r fc belonging to [0,6) C F^Q^ ,6)) n FT" ((a, £7]), (5 > 0, we may set £^(r fe ) = 
(F+)- 1 (r i: ) and £7 (r k ) = (FT")- 1 ^). Moreover, 

1 

"Tfc w s (£ fe (rfc)) - a fc 

This formula combined with monotonicity of u> s yields, 



< < +77+7 ^ . ~^77 < 



a k - w s (b) dr k wf(^(T k )) - a k ' wf(a) - a k 



(r k ) 



< 



w s (£ fe (rk)) - a k 
(5.10) 

for a.e r k . If tt) s (£^) / afe, then it follows that £, k (-) is Lipschitz continuous on [0, e], for some e > 0. 
A similar statement is valid for £7 (•). 

(b) Functions £j (77, 77+1), . . . , (j£, r _i(ri+r-x, r; +r ) are defined as unique solutions to the decou- 
pled system of linear equations (15.71) for any given 77, . . ., 77 +r . This is indeed possible because ak 7^ 
oik+i- The solution £^" depends linearly upon r k , r k+ \. Subsequently, we set £7 +1 := £fe~( r fc> Tk+x), 
k = l,...,l + r-l. □ 

Remark. In the case (a) the derivatives are never zero. They may converge to infinity at t = t(, 



as well as at t = U + t*, if at that time instance w s (^) = a k . 

The lemma above expressed the evolution of the pre-images of facets in terms of r k , i.e. the 
amount of vertical shift of the line l k (-, £~ , w(£,~)). However, in order to render (15.31 ) meaningful, we 
have to figure out the time dependence of r k . At the same time we have to construct fi. We begin with 
an explicit case. 
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Lemma 5.2. Let us suppose that Fk(£Z,£t) is neither of zero-curvature nor interacting and it may 
be degenerate. Then there exist S17, Sl^ £ dJ(ctk) and a unique solution r k : [t*,t* + T max ) — ► R to 
the equation 

~7T = 7+7 x jzr, o r fe (t*) = 0. (5.11) 

They are such that the function 

n(x,t) = "» ~" fc (s - £(t)) + ^ (5.12) 



and to defined by &5.3\) satisfy 

-Q^(s,t) = -fo(s,t) fort £ [t*,U + T max ), s £ (f* (7*(t)),#fa(t)). 

Proof. The non-interaction assumption implies that 

a fc _i = a fc - Aa < Wq s (C) < a k < w^ s (£ + ) < a k+1 = a k + Aa 

or 



Keeping this in mind we set 



Sit = lim -z-(wo,x( x )) = lim Tr(™o,x(s)). (5.13) 

♦(#(*))+ ^ 



X- 



Of course O^", fijj" £ dJ(a k ). We notice that both quantities are well-defined for regular as well as 
degenerate facets. 

Now, we turn our attention to equation (15 -lib , we notice that this equation states that the time 
derivative of r k equals the slope of the straight line passing trough the points (£~ , f2 (£7 ) ) and (£T , SI, (^t ) ) 
. This line provides a section of d J, necessary to construct solutions to (I2.ll ). 

The numerator of (15.1 lb is constant and if £7" (•), £!"(•) are Lipschitz continuous and £Z(rk) > 
£Z(t]c) for all the values of r k , then (15.111 ) has a unique solution. If however, £^(0) = £IT(Q), 
then (15.1 11 1 is singular and this equation requires special attention. A similar situation arises when 
w s(£,^) = Oi k - Fortunately, due to a simple structure of (15.1 II ) we may resolve these issues. 

The ODE ( 15.111 ) governing the behavior of a non-interacting facet F k is obtained by taking the 
time derivative of (I5.5I ). 



^M£(t)) " «k#(*)) = J t Tk(t), J t (^k(t)) ~ = Jfkit). (5.14) 



In reality, we do not assume that w is differentiable everywhere, but its one-sided derivatives do exist 

'A: 



at each point. Due to monotonicity of £f the one-sided derivatives suffice in the formula above. 



((F+)- 1 (r fe ) - (F-)-Hr fe ))^ = AO fc . 



By the definition of £v we rewrite (15.111) as follows 

.(In 
dt 

Here, due to the definition of J and (15 - 1 3b . we have 

ASl = ASlk = S1 + -SI" = -. 

2 
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Since £^( r fc) > £ fc (t^) as long as r k / 0, then we deduce that G, the primitive function of 
(F£)~ l (Tk) - (Fj7)~ l (Tk) such that G(0) = 0, is strictly increasing. Thus (15.111 ) takes the form 

| (Gfa)) = Afi 

or G(Tk) = AQt. As a result function is given uniquely by the formula 

n (t) = G~ 1 (AQt) 

andr fe (0) = 0. 

If we now set $7 by formula (15 - 1 2b . then by the convexity of the set dJ(a k ), we conclude that 
t) G dJ(atk). Moreover, for w defined by (15.3I ). the following equality holds by the definition of 

Q, and T k , 

dt [x > t} - dt ~ $®-${t) - dx™ 

fort G [U,U +T max ),x G (^(i)^fc(i))- 

We note that f2, which we so far constructed, belongs to W^QO, 27r)) for each t > i*, if however 
the facet does not degenerate, then Q(-,t*) G W^QO, 27r)) too. 

□ 

We can infer the following observation from Lemma IBTTl and (15- 11b - 

Corollary 5.1. Let us suppose that w s is increasing (resp. decreasing) in a neighborhood of the pre- 
image [£T , £7] of a non-interacting facet. Then, there exists a positive S, such that for t G tjfc + T) : 

(a) if £+ < then i^(n(t)) >5>0 a.e. (resp. |^(r fc (t)) < 5 < a.e.). 

(b) if < £~ , then i^(r k {t)) <5<0 a.e. (resp. (r k (t)) >5>0 a.e.). 

Proof. The chain formula yields ^£, k = a.e. In the case (a), by the geometry of the problem, 

we deduce that ^ > (see (IBTBT) ) as well as ^ > (see (15.1 lb ). Moreover, formulas (l5T6b and 
(15.111 ) imply that none of the factors may vanish, in fact they are separated from zero. 

The remaining cases are handled in the same way. □ 



We shall state a result corresponding to Lemma 15.21 for a set of interacting facets. It will be 
somewhat more tedious. 

Lemma 5.3. Let us suppose that non-degenerate facets Fi, Fi +r , r > interact, while < 
and £jt r < £f +r+1 . Then, there exist continuous functions : [t*,i* + T) — ► R, k = I, . . . , 
I + r, such that they are locally Lipschitz continuous on (£*, t* + T) satisfying 15.1% below and there 
are C 1 functions r k : [U, U + T) ->■ R, k = I, . . ., I + r, andVL{-,t) G W{(^ (t),^ +r (t)). They are 
all such that w defined by &5.3\) satisfies 



■^W) = fljM) fort G [U,U +T max ), s G (^r(*)>^+r(*))- ( 5 - 15 ) 

Remark. The above Lemma includes the case when the set S\E(wo, s ) consists of a single component. 

Proof. By our assumption the pairs of facets F/_i, F\ and Fi +r , Fi +r+ i do not interact. Thus, the 
evolution of the end points ^ and is determined as for a single non-interacting facet. This remain 
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applicable, unless E(w s (-, ti)) = [0, 2tt). We proceed as in Lemma l5T2l but we have to determine 
Ti + i, i = 1, . . . , r and SI simultaneously. We keep in mind that £r i = £,J +i _ v i = 1, • ■ ■ , r. In order 
to obtain their time evolution, we differentiate (15.71) with respect to time. This yields, 

"fcCfc + h = afe+iCfc+i + f k+i- (5-16) 

The equation for should be similar to (15.111 ). if so we have to select fijlL, i = 1, . . . , r. We define 
and Q7^_ r as in ( 15.131 ), i.e. 

0+ = lim •7-(too,i(a;)), O 

x-»«i r (t))+ ^ 



lim ^-(w 0:X (x)). 

*-($;(*))- dip 



(5.17) 



We have to define the remaining O^'s while keeping in mind f27f = O fc+1 . By the properties of deriva- 



tive %£(<p) and the subdifferential dJ{a.k) the number is one endpoint of the interval dJ(ak), thus 
we inductively define as follows, 



if the facet Ffc has zero curvature 
the other endpoint of the interval dJ(ak), otherwise. 



We have to check that Or\_ r defined in this way agrees with (15. 17b ?- We prove this by induction with 
respect to r, the number of interacting facets. If r = 1, then the claim follows from the preceding 
considerations. Let us suppose validity of the claim for some r > 1, we will show it for r + 1. Let 
us suppose that u>o corresponds to a group of r + 1 interacting facets satisfying the assumptions of 
the Lemma. We consider such a mollification Wq of wq in a neighborhood of £j +r = £^ r+1 that 
Wq = wq for x satisfying \x — £t r | > e and u>q is smooth. Moreover, we require that wo tS and Wq s 
are simultaneously increasing or decreasing. Thus the facets corresponding to Wq are F\, . . . , F/ +r , 
^i+r+l- We notice that facet i^+ r +i is of zero curvature iff facet Fi +r+ i is of zero curvature, i = 0, 1. 
Moreover, facets -F}+ r , -F}+ r +i do not interact. By the inductive assumption £1^+1 = is equal 
to At the same time = 0~ j , 1 * s determined from fiL-/ , 1 = £1^+1+1 an d ^o,s as in 

Lemma I5T21 The two ways of course coincide, due to formulae (15 - 1 3b - Our claim follows. 
We now write equations for k = I, . . . , I + r, they are as (15-1 lb . 



dt 



07 







for k = I, . . . ,1 + r. 



(5.18) 



Since we do not admit degenerate facets, these equations are not singular. We combine them with 
( 15.131 ) and after writing rj = . . . , Ci. r _i)» we arrive at 



Ar) = B(ri), 



(5.19) 



where 



A 



- n 



o 



-oci+i 

(Xl+l 



-Otl+2 



n + — n~ 






k = l + l,...,l + r-2, 
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d iX ^+r-l , ^l+r ^l+r 



B{V)l+r-l = ai +r — il + + + 



Under our assumptions, there is a separate equation for Jj£t r i.e. d5.11| >. Due to the assumption of 
absence of degenerate interacting facets, this system is uniquely solvable on [i*, t* + T). 
We have to define Q, it will be a continuous piece-wise linear function, 

n(*,t) = " (x - er+i(i)) + <v ( 5 - 2 °) 

Moreover, 10 and satisfy (15 -15b - □ 
We claim in theorem 15. lt hat the number of facets decreases in time. The result below explains 
that certain phenomena are forbidden. Namely, no facet with non-zero curvature may degenerate. 

Proposition 5.3. In any group of interacting facets Ff,, k = I, . . . , I + r, r > only a facet with zero 
curvature may degenerate. 



Proof. Let us suppose that Fi, . . ., Fi +r , r > is a maximal group of interacting facets with non- 
zero curvature. For the sake of definiteness, we will proceed while assuming that w s is increasing on 

Step 1. Let us observe that for a facet Fk to disappear, it is necessary, (but not sufficient) that one of 
neighboring facets moves upward faster than i.e. either Vk+\ = rfr ( ^ +1 > = Vk or V^-i = 
> = Vfc. Indeed, the position of Ft is defined by the intersection of the lines containing 
Ffc, Ff. + x moved vertically by and respectively by Tfc +1 and the intersection of lines containing Fk, 
Fk-i moved vertically by Tfc and respectively by Tf--\- Thus, if the lines containing F^+i and F^-i 
are moved up so much that their intersection is above the line containing Ff. moved vertically by T&, 
then facet Fk is going to disappear. This situation may occur only if Vk+i > Vk or Vk-i > Vk- 
Step 2. Let us suppose that facets Ff., F^-x interact, hence by (15.191 ) 

ak£,k - oik-iik-i = Tk-i - f k . (5.21) 

By the monotonicity assumption on w s we notice that fk-i and fk are positive. If the length of Fk, 
which is equal to — stays bounded on [t*,t* + T) while the length of F k -x vanishes at 

t = U + T, then in a neighborhood of t* + T we have ik-x — fk < 0. Thus, by (15.211 ) we can see that 

afc-i(4 - £fc-i) + (afc - afc-i)4 < 

and by ( 15.181 ) the left-hand-side (LHS) converges to — oo when t tends to t* + T. Since — 
must be bounded from above, we deduce that < for i close to t* + T. 

Sfe/? 3. Since always < and > (unless H(</>) = [0, 27r)), we conclude that not all 

of the facets vanish simultaneously at t = t* + T. As a result we may assume the length of 
is greater than d > on [ij, T). Thus, we conclude by step 1, that for £ close to + T we have 
V/ > V/_i. By induction we obtain that 

V k+1 >Vk, k = l,...,j + r-l. (5.22) 

We notice that we have the following possibilities for facet Fi +r : (a) there is an adjacent zero- 
curvature facet F; +r+1 ; (b) £^ r+1 is defined as (i ? ^ r+1 ) _1 (T; +r+ i) (see the proof of Lemma l5.il ). In 
case (a) we can see that Ti +r+ \ = while in (b) r; +r+ i > 0. 
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The condition (15.221 ) combined with (15.181) implies that 



%l €l > ■ ■ ■ > C+z+i C+i+i- 
Hence, the endpoints of Ff., k = I, . . . , j + r — 1 converge to a common limit p. But by step 2 

This is a contradiction, our claim follows. □ 
This observation shows that the initial time to = is special. If the data are poor from the view- 
point of dynamics, but still acceptable, then they get immediately regularized. That is all non-zero 
curvature degenerate facet become regular. 

We are now ready for the proof of the main result. 

Proof of Theorem \5.1\ Part A. We start with data free from degenerate interacting facets. We set 
to = 0, we have to define time instance ti, i = 1, . . . , M postulated by the theorem. We shall proceed 
iteratively. 

It follows from Proposition 15.31 that degenerate, non-zero curvature facets are possible only at 
t = 0, i.e. at the initial time instance. 

Let us suppose that (ti), (ti)] is a connected component of Ei(w s (ti)). We have six possi- 
bilities for F k = F k (£(ti),t£{U)y. 

(a) Fk is regular, does not have zero curvature, is non-interacting; 

(b) Fk is regular, does not have zero curvature, is interacting; 

(c) Fk is regular, has zero curvature, is non-interacting; 

(d) Fj, is regular, has zero curvature, is interacting; 

(e) Fk is degenerate, does not have zero curvature, is non-interacting; 

(f) Fk is degenerate, has zero curvature, is non-interacting. 

Cases (a) and (e) are solved in Lemma l5T2l where corresponding £r are constructed. 
The construction of g£ corresponding to (b), (d) is performed in Lemma [531 We stress that in all 
these cases Tf., is given by (15.111) . 

The definition of is simple if (c) or (f) holds, we just set 

6"(*) = £fc. = r k (t) = 0. (5.23) 

We have to define 11 By the very definition of zero-curvature facets the intersection dJ(£^ + e) n 
dJ(^ — e) is a singleton {a} for any positive e < min{^ fc+1 — — £ k _i}- Moreover, a G A, 

hence we set 

n(x,t) = a, for x G (5.24) 

Thus, we have specified evolution of for every configuration. In all these cases the functions 
k = 1, . . . , Ni are defined on maximal intervals [ij, ti + Tj^}. The numbers are defined as 
follows. 

In (a) and (e) the positive number (resp. T^) is such that < (resp. < 

Cfe (*)) for t <U + T k + (resp. t < ti + TjT), while equality occurs at t = Tjt (resp. t = T k ~), i.e. the 
facet begins to interact with its neig hbor. By Corollary are finite. 

If a group of interacting facets Fi, ...Fi +r does not contain any zero-curvature facet, then by 
Proposition 15.31 it may not vanish and its maximal existence time is defined as in (a) for £/ +r . Thus, 
at T+_ l the group begins to interact with another facet. On the other hand, if this group of interacting 
facets Fi, ... Fi +r contains a zero-curvature facet, say F p , then is defined as the extinction time of 
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F p , i.e. £ p (t) < £ + (t) for te[U,U + T+), while f~ (U + T+) = + T+). Thus, the number of 
facets drops by one. 

Cases (c) and (f) do not contribute to the definition of t i+ i, because ( 15.231 ) is valid for all t>ti. 

We have to define also U(x,t). An attempt to do so reveals another difficulty related to con- 
struction of £r starting from t = 0. Let us consider two interacting facets F).(£7 , Fk'i^k'^k')' 
where 

CH Z K, S ) ; [4",4'l cS r(%). (5-25) 
It is obvious that for any s E [£7 , K"] and s' E the intersection 

5Jh, s M)naJK s (s')) 

is non-empty if and only if \l — r| = 1. If the above intersection is non-empty, we can construct the 
desired Q(x, t). On the other hand, if this intersection is void, then we have no chance to construct a 
Wl section of dJ(w s ). 

Let us suppose then that (15.251 ) holds and \l — r\ = p + 1, p > 0. Let us suppose for simplicity 
that I < r. Thus, a single point £ is a connected component of Hj(ii?o )S ), j = U + 1, . . . ,r, i- e - 

e=$7=#, j = U + l,...,r. 

In other words, we have a number of degenerate, interacting facets at £. The system of ODE's ( 15.191 ) 
is singular. The problem of evolution of interacting degenerate facets shall be dealt with below in Part 
B of the proof. It occurs only at t = 0. 

Finally, we check that w(x, t) and Q(x, t) fulfill the conditions postulated in the definition of the 
weak solution. They satisfy the equation 

w t (x,t) = Q x (x,t) (5.26) 
and the initial and boundary conditions are satisfied. Integral identity in Definition 12. 1 I follows. 

Part B. After finishing part A, i.e. the case of data satisfying (15 -4b . we consider the interaction 
of facets during creation, i.e. (15.41 ) is no longer valid. We have two cases to consider: 

(g) Fk is degenerate, with nonzero curvature and interacting; 

(h) Ffc is degenerate, with zero curvature and interacting. 

We begin with (g). Let us suppose that u>o violates (15.41 ) at some £. Thus, we are dealing with the 
situation when one sided derivatives of wq differ at £, i.e., 

w 0,s(0 <»k< Wo,s& 

for some E A. It may as well happen that the reverse inequalities occur, however for the sake of 
defmiteness we shall stick to the above choice. 

We shall construct two functions w e , w e such that their derivatives belong to J-R, w € (x) < 
wo(x) < w e (x) and 

\w e (x) — wo(x)\, \w e (x) — wq(x)\ < e. (5.27) 

We set 

w e (x) = max{wo(x),l k (x,£,w (Q + 5)}, 
where 5 > is so chosen to guarantee (15.27I ). We also define 

w t (x) = max-OoOr) - e, l k (x, £, itf (£) + 6)}, 
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where 5 G (0, e) is arbitrary. Of course, (15.271 ) holds. 

If the newly constructed w e and w e do not satisfy (15.41 ). we repeat the above process until they 
do. Subsequently, we apply the results of Part A to w e and w e . We deduce from that existence of 
interacting facets at £. By the comparison principle, the non-zero interacting facets exist for t > to- 

Finally, we study (h). We notice that at such an instance F k cannot interact with two neighboring 
facets, because this would mean that F k ^\ and F k+ i lay on the same line, that is, F k is their common 
end point. Thus, the three facets F k ^\, Ft and F k+ \ form a single facet F k with the pre-image 
[£fc_i> £fc+lJ - ® n tne otrier hand it may happen that F k is a degenerate, zero curvature facet interacting 
with just one neighbor, say F k +i- Since Ff. is degenerate, i.e., £t = £r =: due to its interaction 
with F k+1 we have = £ fc . Moreover, w^ s (£ k ) = a k ^ x and u>j£ s (ffc) = a k where a fc _i, a fc G 
.4 and we may assume that a k -i < a k , (the other case is handled similarly) and w s (£t_ 1 ,£i t ) C 
(afc_2, afc-i). A similar situation occurs when interacts with i*V_i. 

In order to determine the evolution of the system we have to take into account if F k+ \ has zero- 
curvature or not. In the former case f k +i = 0, hence we set r k = 0. In latter case we have f k+ i > 
(it may not occur fk+i < 0). Thus, F k disappears instantly. As a result, we agree to disregard F k and 
diminish iVo by 1 . 

Part C. We have to deal with the points outside of E(w s (-,t)) = \J^ X , By the definition 
of E(w s (-, t)), its complement is open 

Ni 

[0,27r]\HK(-,t)) = |J(e+ e fc+1 ). 

1=1 

where , ^ i+1 ) should be understood as , 2ir] U [0, £± ), (with the understanding that < £ fe < 
27r, fc = 1, . . . , JVj). Using again the definition of S, we come to the conclusion that, if x belongs to 
any of the intervals £;T+i)> then either t) exists or wf(x, t) ^ w~ (x, t). In either case, the 
set dw(x, t) (see (14.51 )) does not intersect A. Since dw(x, t) is an interval, we deduce that there exists 
a k e A such that 

dw(x,t) C (a k ,a k+1 ). (5.28) 

We have to make sure that the choice of a k , in the formula above, depends only on the interval 
b ut it is independent from a specific point x G (Cfe )Cfe+i)- Indeed, by the definition of 
the J-R class dw = M — f or <9w; = / — M , where / is a continuous increasing function and M a 
maximal monotone operator. Thus, the images /(^T, Cfe+i) an d j are connected intervals, 

so is the image £jfc+i)> which is disjoint from A. Our claim follows. 

As a result, our definition of w[x, t) for x G" E(w s (-,t)) is as follows, 

w(x,t) = w(x,t k ) and n(x,t) = -j^(w s (y,t k )) for a; G Ci+i)- 

where y G (i^ 1- , i^+i) * s anv differentiability point of w(-, t k ). 

Part D. We have to define t k+ i. We do this inductively. Once t k is given, we set 

tfc+l = tfc + min{min^ + ,min^"}. 

Thus at t k+ \ two facets begin to interact, due to the shrinkage of to a point or due to the 

disappearance of a facet. By Proposition 15.31 we know that only zero-curvature facets may disappear. 
We set 

N i+1 = Ni-m, 
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where m is the number of removed degenerate, interacting, zero-curvature facets at t = U+i. 

The last thing to show is the estimate ||io s (-, £)||j-r < ||w s (-, s))||j-r, whenever t > s. By 
the construction above, the number of connected components of E(w s (-,t)) drops at time instances 
tk, k = 1,..., Mn, hence K(w s (-,t)) < K(w s (-, s)), whenever s < t. It remains to show that 
\\ w s('i t)\\TV(S) — \\ w s('i 3 )\\tv(S)> where we denoted by the total variation of function / 

over set E. 

We first consider the case t > s such that E(w s (-,t)) / S, we know that we always have 
E(w 8 (-,t)) D E(w s (-,s)) for s < t. By the general properties of the total variation, we notice 
that 

\\ w s(-,t)\\ TV (s) = \\Ws(;t)\\TV(S{t)) + \\Ws(;t)\\TV(S\B(t)), 
where we wrote S(cr) for E(w s (-, a). Now, by the definition of w(x, t), we notice that 

lk«(->*)llrv(S\s(t)) = \\w s (-,s)\\ TVi s\s(t)) < \\ws(-,s)\\ TV{s \s(s))- 

We turn our attention to ||ws(-,i)llTV(S(t))- O n tne intervals forming E(t) function w s (-,t) is piece- 
wise constant. The jumps occur at the endpoint of these intervals. They are no bigger and no more 
numerous than the jumps of w s (-, s). Thus our claim follows in the considered case of t. In fact, the 
case of t such that H(t) = S is not much different. Finally, we can see that w s is a difference of two 
monotone functions and one of them is continuous, the other one a maximal monotone operator. 
Our theorem is proved. □ 
We close this subsection with a formula, which might be called "morphing a circle into a square". 
Example. Let us suppose that 4>o(s) = s or wq(s) = \s 2 . Due to the high symmetry of the problem, 
it is sufficient to consider just formation of one facet. Then, w(x, t), the unique solution to (12.11) . is 
given by the formula, 



w(x, t) 



2 



a€[0,£-(t)]U[tf(t),§] 



Here, £f = \ ± <JW X and n = (^-t)^- Let us note that at T\ = ir 2 /2 6 we have 0+ - J2~ = 
^ + — so for later times f = 1. 

We can make this observation more general. 

Proposition 5.4. Let us suppose that the assumptions of Theorem 15. J I are satisfied. Then, there 
exist Tf a , such that if t > Tf a , then w(-,t) is fully faceted, i.e. w s (-,t) is piece-wise linear. More 
precisely, forZ(w s {-,t)) C [0, 2vr) fort < T fa and~(w s (-,t)) = [0,2vr) fort > T fa . 

Proof. Let us consider wq. It is fully faceted or not. If not, then by the proof of Theorem 15. II we 
deduce that after at some tj we have S(tj ) = [0, 2ir) and our claim follows. □ 

5.3 Convexification 

We show that after some depending upon the initial data, the solution becomes such that w s = tp is 
monotone decreasing or increasing. We shall call this process by convexification. 

Proposition 5.5. Let us suppose that the assumptions of Theorem I5.il are satisfied. Then, there 
exist T cx , such that if t > T cx , then w s (-,t) is monotone, while this is not true for t < T cx . 
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Proof. If WQ tS is monotone, then we are done. Otherwise, let us suppose that tj is the largest time such 
that at tj a zero curvature facet disappears. Since the zero-curvature facets cannot persist because their 
endpoints necessarily move, it follows that T cx = tj has the desired properties. □ 

Remark. All possibilities can be realized T cx > Tf a as well as T cx < Tf a . 



5.4 Asymptotic behavior of facets 

Here, we consider the last stage of evolution, when t >tu and Nm = 4. In this case, it is sufficient to 
specify only £f, k = 1,2, 3, 4. Furthermore, the system for interacting facets, (15.191) takes the form, 



(5.29) 



We notice that the stationary points of ( 15.291 ) are such that f\ = ... = T4. This occurs if and only 
if $7^ — = ik — k = 1,2,3,4, where by £0 we understand £4. Moreover, due to our 

assumptions on J we have — = AQ, k = 1, 2, 3, 4. 

Additionally, system ( 15.291 ) possesses a Liapunov functional. Namely, let us write 



- 


026 


= Tl - 


T"2 


026 - 


«3?3 


= 72 ~ 


T"3 


036 - 


«4?4 


= T 3 - 


T"4 


a 4 £ 4 - 




= f 4 - 


Tl 




= 


, fc = 


1,2,3,4. 



F(0 = ]Tln(£ fc -6- 



lj 1 



fc=l 

with the understanding of £0 as above. By direct calculation, we check that 

This derivative vanishes if and only if £ is the only equilibrium point. Thus, we have a complete 
picture of the asymptotic behavior of A. 



Theorem 5.2. Let us assume that ipo £ J-R and w is the corresponding unique solution to i2.1\) . 
Then, there exists T\, max{T cx , Tj a } < T\ < 00 with the following property: 

(a) IfTi < 00, then £ t (t) = -^f + f I + a, for some a > 0, 1 = 0, . . . , 3, and t > Ti, in other words, 
w is the minimal solution for t > T%; 

(b) If Ti = 00, then lim^oo f,(t) = + f I + a, I = 0, . . . , 3 for some a > 0. 

□ 



5.5 Proof of Theorem |23] 

In the course of proof of Theorem l5.1[ we exhibited a quite explicit construction of the weak solution 
with such initial data that ipo £ J-R. Now, we have to show that is has all the postulated properties 
of the almost classical solution. We have already noticed that w s = A s + s belongs to the J-R class, 
furthermore \\w s (-, t)||j-R < \\w s (-, 0) || j-r- The key point, however, is to realize that 

Q = dJodw, (5.30) 
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where dw is the multivalued map whose section is w s . We defined dw in (14.51) . Checking that (15.301) 
indeed holds requires recalling the steps of construction of S7, we will do this below. Finally, after we 
set N = {0, tx, . . . , tjvf} we see that 

d 

A t = —dJo(A s + s ), 

holds for all t £ (0, +oo) \ N in the L 1 sense, more precisely it holds pointwise except x G 
[0, 2vr) : i = 1, . . . , iVfc}- Indeed, the definitions (f57TTb . (f57T8T> . (15331 of Tfc(i) were such 

that JjTfcO) = f^. Moreover, ^ = ^r fc (t), see LemmaEU Lemma|53]and eq. K26\> . We recall 
that by definition functions Tfc(-) are continuous on [i,, and differentiable in (ij, tj+i). Moreover, 
the right derivative of Tfc(i) is well-defined for all t, except possibly t = to. Hence, ^ is defined 
everywhere, except the points ij, % = 0, . . . , M, but the right time derivative is defined for all 
t > 0. 

We will check below that f2, constructed in the course of proof of Theorem 15- 1 L coincides with 
dJodw, — see (15 - 1 2b . (15.201 ). (I5.24I ). where w s (s, t) = A s (s, t) + s. In order to see that we examine 
the steps of the construction of f2 and compare it with the definition of the composition o. Let us fix 
t € (tk,tk+i), at the end we will consider t = tk+i, then we compose dw(-, t) : [0, 2tt] — ► [a, b] with 
d J : E — > M. We have to identify the sets V s , Vf and P r appearing in the Definition |2.3i For our 
choice of t we have 

^) = {«e [0,2tt] :w+(s,t)^w-(s,t)}. 
In particular, P s (t) contains all points £f(t), i = 1, . . . , -/V^. We can see that 

i=i 

i.e., it is the sum of interiors of intervals contained in E(w s (-,t)). Finally, by the definition 

D r (t) = [0,27r]\(D a (t)UP/(t)). 

We shall consider these cases separately. 

1° case V r . If s 6 V r (t), then u; is differentiable at s and w s (s,t) Thus, by ( 12.41) 
dJ5dw(s,t) = ^(w s (s,t)). We notice that P r (t) C [0,2vr] \ E(w s (-,t)), hence by Part C of 
the proof of Theorem 15. II we immediately see that dJodw(s, t) equals Q(s, t) on V r . 

2° case Vf. By its definition T>At) is the sum of interiors of pre-images of facets, as noticed 
above. Moreover, on each interval (£r (t),^(t)), the set dw(x, t) is a singleton equal to {ctk} C .A. 
Then, the cases of the Definition 12.31 see formulas (I2.5I )- (I2.8I ) have their counterparts in the formulas 
Km . KM and (EH- 

3° case V s . We notice that, if t > 0, then the set 3(iu s (-,i)) has no component, which is a 
singleton. Thus, if s G T> s (t), then the set dw(s,t) does not intersect A. As a result, formula ( 12.91 ) 
for the composition yields a singleton, because on the RHS of (12.91 ) the limit of constant functions are 
taken. This in agreement with the discussion of Part C. 

Finally we have to deal with the case t = ifc+i. On one hand tk+i) is defined by the left time 
continuity of Q, on the other hand we have to check that f2 = dJodw. 

By the very definition of tk+i (see Part D of the proof of Theorem 15.11) . at this time instant a 
zero-curvature curvature facet disappears or two facets begin to interact or merge, i.e., 

lim t+(t) = a= lim £r+iW- 
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We have then two possibilities, either a = £J(t k+1 ) = £ j+1 (t k+1 ) or a G (£fc+i),£t(*fc+i)) 
where this interval is a connected component of E(w s (-, tfc+i))- Once we realize this, it is clear that 
t k +i) = dJodw(-,t k+1 ). □ 



6 Appendix 

6.1 Motivation of equation dl.2D 

Here, we consider closed curves, we view them as graphs over a smooth, convex reference closed 
curve M.. We do not make here any attempt to consider non-smooth reference curves, which is 
reasonable because this would add up difficulties while not giving advantages. 

Let us suppose that xq(s) is an arc-length parameterization of M. and e t (s), e n (s) are unit tangent 
and normal vectors, respectively, such that (e t (s) , e t (s)) is positively oriented. Then all points in a 
neighborhood of M. can be uniquely written as x = xq(s) + e n A, as a result we can parameterize our 
curve r(t) as 

x(s, t) = x (s) + e n (s)A(s, t). 

Since M. is convex we may write e n uniquely as e n (ip(s)) = (cosip(s),smip(s)), where tp is the 
measure of the angle between the x\ axis and e n . Moreover, 

^e n (ip(s)) = -e t (ip(s))^ = -Ket(<p(s)). 

We note 

— (s,t) = e t (l - kA) + e„A s , 

because |xo(s)| = 1- With this formula at hand, we can write the expression for the tangent and 
normal to T(t), they are r = pp(e t (l - kA) + e n A s ), n = ^(-A s e t + (1 - nA)e n ), where 
W 2 = (1 - kA) 2 + A 2 S . Hence, the LHS of 0> takes the form 

(3V = l3^-n=±-(l-KA)A t . 
at W 

The RHS of dl.ll ) is div5V^7(^)|^ =n . In our paper |MRy[ , we have shown that it is equal to 

k = is ih Ie{lf) 

We defined I$(<p) as follows, I&(<p) = j(n(<p)) + ff dip ft j(ia.(t))dt. We noted that this function 
is convex iff the stored energy function 7 is convex. However, in general 1$ does not enjoy higher 
regularity properties. It is not differentiable at angles corresponding to the normals to the Wulff shape. 
Finally, equation (11.11 ) takes the form 

0n-e n A t = ^(J^Ig(a)\, (6.1) 

where a is the measure of the angle between the x\ axis and n. 

One may study evolution of convex curves defined by their angle parameterization. We notice 
a = <p+ip, where tp is the measure of the angle between r and e t . We notice that r -e n = sin ip = 
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r-e n = cosip = -^p. Thus, we see that ip = Arg(r-e n +iT-e n ), in fact we have ij) = arctan 
Thus, (11.11) takes the form 

d ( Q T ( /A, 
pn • e n A t = — -rr-le f + arctan 



ds \d4> ° V V 1 - A 

This equation is rather involved, we prefer to simplify it by dropping the terms which at this stage we 
deem not important, thus we come to (1 1 -2b - 



6.2 Other choices of function J 

We may also consider any properly chosen piecewise linear, convex J, 

N 



Ji{<£) = S ^b i \ L P ~ a i\- (6-2) 



8=1 

We require that N > 4, b{ > and cto < a\ < ... < < ao + 2ir, we will write S = [ckq, ao + 2-7r). 
In order to stick to geometrically relevant data, we also impose the condition that ^2u=i = n > which 
guarantees that dl(S) is an interval of length 2tt. In addition, we assume that the following function 
yields an angle parameterization of closed curve, which encompasses a convex region. Namely, we 
set 

3 n 

^ = Y. bi - Yl 6i > 3 = 0,...,N, (6.3) 

i=i *=i+i 

with the convention that the summation over an empty set of parameters yields zero. Then we define 
$ : [ao, «o + 2vr) — > R by the formula 

N 

8=0 

(with the convention «at+i = ao + 2vr) is an angle parameterization of closed curve. We notice that 
our assumptions imply that Qq + 2ir = 0,^. 

The analysis of behavior of solutions presented in Section [5] is valid also for J given by dl.4| ) and 
Jl, however the actual calculations for J\ are more lengthy. In addition we may show existence of 
weak solution for a general, piecewise smooth, convex J, but in this case we cannot offer detailed 
analysis of solutions, yet. 
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